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Abstract 

We develop a theory of local asymptotic normality in a quantum regime based on a novel 
quantum analogue of the log-likelihood ratio. This formulation is applicable to any quantum sta- 
tistical model satisfying a mild smoothness condition. As an application, we prove the asymptotic 
achievability of the Holevo bound for the local shift parameter. 

1 Introduction 

Given a (classical) statistical model S = {pe ; 6 G Q} on a. probability space (17, J^, ji) indexed by a 
parameter 6 that ranges over an open subset Q of M'', let us introduce a local parameter h := ^/n{9—9o) 
around a fixed 9o G Q. If the parametrization 6 i-^ pg is sufficiently smooth, it is known that the 
statistical properties of the model jp^'^^y^ ; /i e ]R''| is similar to that of the Gaussian shift model 
{N{h, JgJ^) ; h € R''} for large n, where " is the nth i.i.d. extension of pg, and Jg^ is the Fisher 
information matrix of the model pg at Oq- This property is called the local asymptotic normality of 
the model S pU] . 

More generally, a sequence of statistical models on (fi'"^ J^' "),^(")) is called 

locally asymptotically normal (LAN) sd 9q €z Q if there exist el d x d positive matrix J and random 
vectors A(") = (aJ"\ . . . , A|;"') such that A^") ^ N{0, J) and 

(») 

log ^^n^ = - l^^^^^ + 

Peo 

for all heR'^. Here the arrow -^^ stands for the convergence in distribution under pi"-', , , ^, the re- 

mainder term Op^^^ (1) converges in probability to zero under and Einstein's summation convention 
is used. The above expansion is similar in form to the log-likelihood ratio of the Gaussian shift model: 
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This is the underlying mechanism behind the statistical similarities between models ^Pg^^f^/^ ; ^ £ M.'^"^ 
and {N{h,J-^)] h&R"^}. 

In order to put the similarities to practical use, one needs some mathematical devices. In general, 
a statistical theory comprises two parts. One is to prove the existence of a statistic that possesses 
a certain desired property (direct part), and the other is to prove the non-existence of a statistic 
that exceeds that property (converse part). In the problem of asymptotic efficiency, for example, the 
converse part, the impossibility to do asymptotically better than the best which can be done in the 
limit situation, is ensured by the following proposition, which is usually referred to as "Le Cam's third 
lemma" [20] . 

Proposition 1. Suppose is LAN at 9o G Q, with A^"^ and J being as above, 

and let X*-"-* = (X^"\ . . . , X^"^) be a sequence of random vectors. If the joint distribution of X^""^ and 
converges to a Gaussian distribution, in that 




then ^ N{Th, E) for all /i e M''. 

Now, it appears from this lemma that it already tells us something about the direct problem. In 
fact, by putting := J^'Li [J^'^Y'' A^"\ we have 




so that X^"^ N{h, J~^) follows from Proposition [T] This proves the existence of an asymptotically 
efficient estimator for h. In the real world however, we do not know Oq (obviously!). Thus the existence 
of an asymptotically optimal estimator for h does not translate into the existence of an asymptotically 
optimal estimator of 6. In fact, the usual way that Le Cam's third lemma is used in the subsequent 
analysis is in order to prove the so-called representation theorem, [20l Theorem 7.10]. This theorem 
can be used to tell us in several precise mathematical senses that no estimator can asymptotically do 
better than what can be achieved in the limiting Gaussian model. 

For instance. Van der Vaart's version of the representation theorem leads to the asymptotic minimax 
theorem, telling us that the worst behaviour of an estimator as 9 varies in a shrinking (1 over root 
n) neighbourhood of Oq cannot improve on what we expect from the limiting problem. This theorem 
applies to all possible estimators, but only discusses their worst behaviour in a neighbourhood of 9. 
Another option is to use the representation theorem to derive the convolution theorem, which tells us 
that regular estimators (estimators whose asymptotic behaviour in a small neighbourhood of 9 is more 
or less stable as the parameter varies) have a limiting distribution which in a very strong sense is more 
disperse than the optimal limiting distribution which we expect from the limiting statistical problem. 

This paper addresses a quantum extension of LAN (abbreviated as QLAN). As in the classical 
statistics, one of the important subjects of QLAN is to show the existence of an estimator (direct part) 
that enjoys certain desired properties. Some earlier works of asymptotic quantum parameter estimation 
theory revealed the asymptotic achievability of the Holevo bound, a quantum extension of the Cramer- 
Rao type bound (cf.. Appendices El [B]). Using a group representation theoretical method, Hayashi and 
Matsumoto [8j showed that the Holevo bound for the quantum statistical model 5(C^) = {pe ; 9 G 8} 
comprising the totality of density operators on the Hilbert space "H ~ is asymptotically achievable 
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at a given single point Oq E Q. Following their work, Gut^a and Kahn [51 [TT] developed a theory 
of strong QLAN, and proved that the Holevo bound is asymptotically uniformly achievable around a 
given 9o G Q for the quantum statistical model iS(C^) = {pg ; 9 G 0} comprising the totality of density 
operators on the finite dimensional Hilbert space H ~ C^. They proved that Ipf""^ ; /i € M''! and a 
quantum Gaussian shift model can be translated by quantum channels to each other asymptotically. 
Although their result is powerful, their QLAN has several drawbacks. First of all, their method works 
only for i.i.d. extension of the totality 3(1-1) of the quantum states on the Hilbert space H, and is 
not applicable to generic submodels of 5(1-1). Moreover, it makes use of a special parametrization 9 
of S(l-l), in which the change of eigenvalues and eigenvectors are treated as essential. Furthermore, it 
does not work if the reference state pg^ has a multiplicity of eigenvalues. Hayashi and Matsumoto's 
formulation j8] also suffers from the same problems. 

The purpose of the present paper is to develop a theory of (weak) QLAN based on a new quantum 
extension of the log-likelihood ratio. This formulation is applicable to any quantum statistical model 
satisfying a mild smoothness condition, and is free from artificial setups such as the use of a special 
coordinate system and/or non-degeneracy of eigenvalues of the reference state at which QLAN works. 
We also prove asymptotic achievability of the Holevo bound for the local shift parameter h that belong 
to a dense subset of . 

This paper is organized as follows. The main results are summarized in Section[21 We first introduce 
a novel type of quantum log-likelihood ratio, and define a quantum extension of local asymptotic 
normality in a quite analogous way to the classical LAN. We then explore some basic properties of 
QLAN, including a sufficient condition for an i.i.d. model to be QLAN, and a quantum extension 
of Le Cam's third lemma. Proofs of those results are provided in Section |3l Section |4] is devoted 
to application of QLAN, including the asymptotic achievability of the Holevo bound and asymptotic 
estimation theory for some typical qubit models. For the reader's convenience, a brief account of 
quantum estimation theory are presented in appendices A-D. Those prerequisites are used throughout 
the paper. 

It is also important to notice the limits of this work, which means that there are many open 
problems left to study in the future. In the classical case, the theory of LAN builds, of course, on the 
rich theory of convergence in distribution, as studied in probability theory. In the quantum case, there 
still does not exist a full parallel theory. Some of the most useful lemmas in the classical theory simply 
are not true when translated in the quantum domain. For instance, in the classical case, we know that 
if the sequence of random variables Xn converges in distribution to a random variable X, and at the 
same time the sequence Yn converges in probability to a constant c, then this implies joint convergence 
in distribution of (A„, Yn) to the pair (X, c). The obvious analogue of this in the quantum domain is 
simply untrue. In fact, there is not even a general theory of convergence in distribution at all: there 
is only a theory of convergence in distribution towards quantum Gaussian limits. Unfortunately, even 
in this special case the natural analogue of the just mentioned result simply fails to be true. 

Because of these obstructions we are not at present able to follow the standard route from Le Cam's 
third lemma to the representation theorem, and from there to asymptotic minimax or convolution 
theorems. 

However we believe that the paper presents some notable steps in this direction. Moreover, just as 
with Le Cam's third lemma, one is able to use the lemma to construct what can be conjectured to be 
asymptotically optimal measurement and estimation schemes. We make some more remarks on these 
possibilities later in the paper. 
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2 Main results 



Definition 2 (Quantum log-likclihood ratio). We say a pair of density operators p and cr on a finite 
dimensional Hilbert space T-L arc mutually absolutely continuous, p ~ a in symbols, if there exist a 
Hermitian operator C that satisfies 

Wc shall call such a Hermitian operator £ a quantum log-likelihood ratio. When the reference states p 
and a need to be specified, C shall be denoted by C {(j\p), so that 

We use the convention that C {p\p) — 0. 

Example 3. Faithful states are always mutually absolutely continuous. In fact, given p > and 
cr > 0, they are related as ct = e^'^'-°'l''^pe5-^'-'^l''\ where 



4P<^4py P M = 2 log a/ y/apy/^ s/^ 



C{a\p) ^2\og 



Note that Trpe^-^'-'^'''-' is identical to the fidelity between p and a, and e^-^'^'^l'') is nothing but the 

1/2 

operator geometric mean a^p^^, where A^B := A^l"^ (^A^^^^BA^^^^) A^^^ for positive operators 
A,B. 

Example 4. Pure states p = \tp) and a = |^) (^| are mutually absolutely continuous if and only if 

m) ^ 0. 

Proof. Suppose first that p ^ a. Then 

I (^1^) |2 = Tr per = Tr | V) (VI e^^^'^l") | V) e^^^'^l"' = (V'l e^^^-^l") |^ 
Suppose next that (^1-0) 7^ 0. Then 



> 0. 



R:=I 



\m 



10 iii - w 



is positive definite, and C {cr\p) :~ 21ogi? satisfies 



\0 



showing that p ^ a. Note that Trpes-'-^'^l'') is the fidelity again. 

Given a. d x d real skew-symmetric matrix S, let CCR(S') be the CCR algebra defined by 



□ 



(see [TllTllTniin])- We caU X = {Xi, . . . , X^) the basic canonical observables of CCR (S). A state on 
CCR (S) is characterized by the characteristic function ^^{f/*} := (/)(e^'^^^'^'), where ^ = (C')iLi G I^'^ 
and Einstein's summation convention is used. A state (j) on CCR (S) is called a quantum Gaussian 
state, denoted by ~ A^(/i, J), if the characteristic function takes the form 



4 



where h = (/ii)iLi G and V = (Vij) is a real symmetric matrix such that the Hermitian matrix 
J := V + \^—lS is positive semidefinite. When the canonical observables X need to be specified, we 
also use the notation {X, </)) ~ N{h, J). 

We will discuss relationships between a quantum Gaussian state on a CCR and a state on another 
algebra. In such a case, we need to use the quasi- characteristic function 



\t=l I \i=l ^ ' t=l s=t+l 



(1) 



of a quantum Gaussian state, where (X, (/)) ^ N{h, J) and {^tl^^i is a finite subset of ^TU\ . 

Given a sequence 'H^"-', n e N, of finite dimensional Hilbert spaces, 
and p^"^ be a list of observables and a density operator on each 'H^'^\ We say the sequence (X^"', p^"') 
converges in law to a quantum Gaussian state N{h, J), denoted as (X*^"-*, /O*-"^) ~~> A^(/i, J), if 



hm Trp(") []e^«;^'"' = m^^'""'' 



\t=i / \t=i / 

for any finite subset {£,t}t=i of C'', where (X, 0) ~ N{h,J). Here we do not intend to introduce the 
notion of "quantum convergence in law" in general. We use this notion only for quantum Gaussian 
states in the sense of convergence of quasi-characteristic function. 

The following is a version of the quantum central limit theorem (see |10j , for example) . 

Proposition 5 (Quantum central limit theorem). Let Ai (I < i < d) and p be observables and a state 
on a finite dimensional Hilbert space % such that Tr pAi ~ 0, and let 



n 



/n 

k=l 



Then {X^^\ p®^) N(0,J), where J is the Hermitian matrix whose {i,j)th entry is given by Jij ~ 
TipAjA,. 

For later convenience, we introduce the notion of an "infinitesimal" object relative to the conver- 
gence (X("), p(")) iV(0, J) as follows. Given a fist = {x[''\ X^"^) of observables and a 

state p^"-' on each "H*^"' that satisfy (X("\ p^^"^) ^ N{0, J) ^ {X, (f), we say a sequence i?^"-' of obscrv- 

q 

ables, each being defined on "H*-"', is infinitesimal relative to the convergence (X("),p(")) ^ N{0,J) if 
it satisfies 



(2) 



for any finite subset of {Ct}t=i of and any finite subset {rit}l^i of C. This is equivalent to saying 
that 



and is much stronger a requirement than 

(i?("\p("))->iV(0,0). 

An infinitesimal object relative to (X("),p(")) -w N{0,J) will be denoted as o(X("), p(")). 

9 

The following is in essence a simple extension of Proposition [5l but will turn out to be useful in 
applications. 
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Lemma 6. In addition to assumptions of Proposition\^ let P{n), n Cz N, be a sequence of observables 
on %, and let 



,1 



k=l 



//lini„_^oo P{n) = and lim„^oo V^Tr pP(n) = 0, then o{X p'^") . 

We are now ready to extend the notion of local asymptotic normality to a quantum regime. 

Definition 7 (QLAN). Given a sequence H^"-' of finite dimensional Hilbert spaces, let 5'"' = 
; ^^ e C M''! be a quantum statistical model on "H^"), where p^"-* is a parametric family of 
density operators and <d is an open set. We say 5^"^ is quantum locally asymptotically normal (QLAN) 
at 6*0 e 6 if the following conditions are satisfied: 

(i) for any e 6 and n e N, p^""* is mutually absolutely continuous to 

(ii) there exist a list A'") = . . . , A^"') of observables on each 'H^"-' that satisfies 

(A("),p(;:))^iv(o,j), 

where J is a d x d Hermitian positive semidefinite matrix with Re J > 0, 



(iii) quantum log-likelihood ratio C^l^^ := C (pi" 



An) 



is expanded in /i e R'' as 



= /i^Af) - i( J,,/i'^/»^)/(") + o(A("), p^^)), (3) 
where /'■"•' is the identity operator on "H*^"^. 

It is also possible to extend Le Cam's third lemma (Proposition [1]) to a quantum regime. To this 
end, however, we need a device to handle the infinitesimal residual term in ^ in a more elaborate 
way. 

Definition 8. Let 5^") = ; 6* e 6 C M'^j be as in DefinitionEl and let = {x'f \ X^"^) 

be a list of observables on H'^"^ We say the pair (iS^"-* , X'")) is jointly QLAN at Oq (z Q ii the following 
conditions are satisfied: 

(i) for any 6* € 8 and n e N, p^""* is mutually absolutely continuous to pg"'' , 

(ii) there exist a list A'") = (A^"-*, . . . , A^"^) of observables on each "H^"^ that satisfies 

(t:V?:')T"((o)-C- 

where E and J are Hermitian positive semidefinite matrices of size r x r and d x d, respectively, 
with Re J > 0, and r is a complex matrix of size r x d. 



(iii) quantum log-likelihood ratio Cj^ := C I p 



is expanded in ft, G M'' as 



r)-i(j,,ft''ft^)/(")+of r!"M,p(:^^ (s) 



With Definition ISl we can state a quantum extension of Le Cam's third lemma as follows. 
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Theorem 9. Let S^'''' and X'^'''> be as m Definition\g is jointly QLAN at Oq e 6, then 



for any h gM."^. 

In applications, we often handle i.i.d. extensions. In classical statistics, a sequence of i.i.d. exten- 
sions of a model is LAN if the log- likelihood ratio is twice differentiable [23 . Quite analogously, we 
can prove that a sequence of i.i.d. extensions of a quantum statistical model is QLAN if the quantum 
log-likelihood ratio is twice differentiable. 

Theorem 10. Let {pe ', E Q (Z K'^} be a quantum statistical model on a finite dimensional Hilbert 
space H satisfying pg ~ pe„ for all 9 G Q, where € & is an arbitrarily fixed point. If Ch 
C {pea+h\pea) *'5 differentiable around h = and twice differentiable at h = 0, then {pf"" ; £ Q C R''} 
is QLAN at Oq: that is, pf" ~ pf", and 



n 

/t7 ^- — ^ 



/n 

k=l 



and Jij := Tr pg^Lj Li, with Li being the ith symmetric logarithmic derivative at 6q S Q, satisfy 
conditions (ii) (Hi) in Definition^ 

By combining Theorem 1101 with Theorem |9l we obtain the following. 

Corollary 11. Let {pe ; ^ £ 6 C M"^'} be a quantum statistical model on % satisfying pg ~ pg^ for 
all 9 Cz Q, where 6q Cz Q is an arbitrarily fixed point. Further, let {-Bi}i<i<r be observables on % 
satisfying Tr pg^Bi = for t = l,...,r. If Ch ■— C (pg^+hlpeo) differentiable around h = and 
twice differentiable at h — 0, then the pair ({pf""} , A"^"-*) of i.i.d. extension model {pf"} and the list 
X^"^ = {X^"'^}i<i<r of observables defined by 

' . An ^^—^ 

is jointly QLAN at do, and 



/n 

k=l 



for any h £ W^, where E is the r x r positive semidefinite matrix defined by Eij = Tr pg„BjBi and t 
is the r X d matrix defined by Tij — Tr pg^LjBi with Li being the ith SLD at 6q. 

As in the classical case, Corollary [TT] prompts us to expect that any estimator for a quantum Gaus- 
sian shift model {iV((ReT) /i, S) ; h€ R'^} could be reahzed asymptotically on jp^^"^^/^ ; h € K'^j. 
This program will be partly demonstrated in Section|3]in the form of achievability of the Holevo bound. 
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3 Proofs of main theorems 
3.1 Proof of Lemma [6] 

We shall prove @ for C C and {vt}t=i C C. 

= Trp«" J[ exp I ^ ^ /^('^-i) (CM. + ^t^(^)) 
= Trp«" n (exp ( ^ feM. + ,7tP(n)) 



i-fc) 



/n 



On' 



Trp' 



nexp(^^(Ct'A, + ,7tP(n)) 
Trp|nexp(^^(CM, + r;,P(n)))j 



/ci , hj- 



kt 



kr = 1 in the summand are 

r 



where Z+ = {0, 1,2,...}. The terms corresponding to fci 
because TrpA^ = and TrpP(n) = o(--^). The terms corresponding to fci + • • • + fc,. = 2 are 

Ul+--- + fcr = 2 \t = l / ) 

t=l t=l ,s=t+l 

t=l t=l s=t+l 



o I - 



t=l s=t+l 



In the third line, we used the fact that P{n) — o(l). Let us denote the terms corresponding 
ki + ■■■ + kr > 3 hy 



r„:=Trp<^ E 



fclH hfcr>3 



^-1 



(feiH hA-v) r 
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Then 



r„ < 



< 



E 

fciH |-A,v>3 

1 



(fel+---+fer) 



n\/n 



E 



fclH |-fcr>3 



t=i *■ 



t=i 



< ^ E n rTll^*^'' + '?*^(")l 

^ fciH hfer>3 t=l ^ 

< ^ E n rrlkt^'^ + ^'^Hl 

^ fci,...,tv6Z+ t=l 
1 '' 

= ^ n '^'^P 11^*^* + VtP{n)\\ 



Since lini„_>.oo P{n) = 0, the operators P{n) are uniformly bounded. As a consequence, hm^ — n [t*^ 
0, so that r„ = o Thus we conchide that 



hm Trp®" TTe' 



t=i 



t=l s=t+l 



( r r r \ 

-2E^j^*'^^-E E ^^'J^A 



The last equation is due to ^ with /i = 0. 
3.2 Proof of Theorem M 

Let Xi, . . . , Ai, . . . , Ad be the basic canonical observables of CCR ( Ini ( ^ J ) ) ' ^^'^ ^ 



quantum Gaussian state N 
quantities 







.T* J J 



R 



on that CCR. Assumption ([5]) guarantees that the 
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enjoy = o y\^(n)j ^ Peaj gslcYi h G K'*. Consequently, for a finite subset {£,tYt^i of C*, 



Trp^+vv^ I n 

= Ir(e^-^'p^)e^-r')(ne^«-^ 



\t=i 



e-^^'^'^-Trp^';) I e 



^(„)^i(/.'A<")+<') l^p^.V-TcJxr'^ ^i(.'A<")+<') 



Since is infinitesimal relative to the convergence (|4]), we see from ([2]) that 
hm Tr o^"^ I TT e^«'-^'"' I 



e 2 '''''' -^'j exp 



V Vt=i / / 

r+1 r+1 r+1 \ 

t=0 t=0 s=i+l / 

\ t=l t=l t=l s=t+l / 

= cxp ( ^ (^/^Cj/i-'" (Rer),,^. - ieje^S,.) -EE ^J^^'^^O ' 

\t=l ^ t=l s=t+l I 

~ It. r\ ~ ~ ~ ~ rr- rr- 

where S :~ and (^o, ■ ■ ■ ,^r, G+i) := (—2^'*' ^i' ■ • ■ 1 —"2^'^)' ^^^'^ © ^^'^ at 

the second equation. This is the quasi-characteristic function of Af((Rer) /i, S). 

3.3 Proof of Theorem [TUl 

Since 



we see that 



-cipf" -E^''*'"'^®^(^''K)®^''*""'^- (6) 
fe=i 

This proves " - " for all 6* G 6 and n e N. 

Before proceeding to the proof of (ii) and (iii) in Definition [51 we give some preliminary considera- 
tion. Let the quantum log-likelihood ratio Ch C ipeo+h\peo) be expanded into 
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where Ai {1 < i < d) and Bjj (1 < i,j < d) are Hermitian operators on H. Observe that Ai is the SLD 
in the ith direction. In fact, 



peg+h = exp 



^{h^A, + o{h)) 



peg exp 



-{h^A., + oih)) 



so that 



So + (Peo^i + Aipeg) + o{h), 



dtPOg = - {peoA + Aipog) . 



This observation also shows that Tr pg„Ai = for ah i. On the other hand, 
Tr pe„+h = Tr pgg exp {h'A, + Bijh'h^ + o{h^)) 

= Trpe„ (^I + {h'A, + B,,h'h^) + i {h'A, + B,,h'h^Y + o{h^) 

= 1 + (Trpeo^«) + h'h^Trpg, {b,, + \A^A^ + oih") 

= 1 + h'^Tvpeg I^B,j + \a,A^ + o(/i2). 

Since Trpgo+Zi = 1 for all h, the above equation leads to 

Trpso (^B,, + ]^A,A^=Q. 

Now we prove (ii). Let Jij Tr pggAjAi, and let 

n 

It then follows from the quantum central limit theorem (Proposition [5]) that (^A^^\ pf") N{0, J). 



(7) 



Finally, we prove (iii). It follows from ([6]) that 



k=l 



Let us show that 



is infinitesimal relative to the convergence (A^"', pf^) N{0, J). It is rewritten as 



" r 1 1 

A.— 1 
n 



J j®(n—k) 



k=l 



^h'A, + -B,,h'h^ + o(-) - -^h'A, + (J„7i*/i^) / 



-1) 



fc=i 



n 

^1 ^ 



(n-fc) 



I j-!gi(n~k) 
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where 

Pin) (^i (^B,, + i J,,/) h^h^ + o(i)) 

Note that hm„_>.oo Pin) = 0, and that 



hra y/nTr p0„P{n) = Tv pg„ ( B,^ + -J^I ) ft' ft-' 



n— J-OO 



= Trpe„ (^By + ft'ft^' 

= Trpeo (^B„ + ^A,A,^ ft^ft^' 

= 

because of ([7]). It then follows from Lemma[6]that r!]^'' ~ o(A("' , pf^) for each h G M''. This completes 
the proof. 

3.4 Proof of Corollary [IH 

That pf" - pf^ was proven in the proof of Theorem [TUl Let aJ"\ . . . , A^"^ be as in the proof of 
Theorem 1101 It then follows from the quantum central limit theorem that 

Further, because of Lemma IHl the sequence of observables given in the proof of Theorem [TU] 
is also infinitesimal relative to the convergence ([8|). Now that (p^", X'"^) are jointly QLAN at 6*0, 
the property ^X^"^ p^J^^^^^ iV((Re T)ft, S) is an immediate consequence of Theorem [9l This 
completes the proof. 



4 Applications to quantum statistics 

Quantum Le Cam's third lemma fCoroUarv fTT|) implies convergence of ft e K''! to a 

quantum Gaussian shift model {A^((ReT) ft, S) ; ft G R''}. This fact prompts us to expect that, for 
sufficiently large n, the estimation problem for the parameter ft of p^^j^y^j ^ could be reduced to that 
for the shift parameter ft of a quantum Gaussian shift model Af((RcT) ft, S). The latter problem has 
been well-established to date (see Appendix [BJ. In particular, the best strategy for estimating the 
shift parameter ft is the one that achieves the Holevo bound C/j (iV((Re r) ft, S), G), (see Theorem 
E5)) . Moreover, it can be shown (see Corollary [M)) that the Holevo bound C/i (A^((Re r) ft, S), G) is 
identical to the Holevo bound Cq^ [pe, G) for the model pe at 6*0 This observation strongly suggests the 
existence of a sequence Af ^"^ of estimators for the parameter ft of that asymptotically 

achieves the Holevo bound Cqo {pe, G). 

This section is devoted to materialize this program: we prove that there exist a sequence of 
estimators on ^^^^ asymptotically unbiased and achieves the Holevo bound Ce^ {pe,G) 

for all ft that belong to a dense subset of R'^. Since this result requires only twice differentiability of the 
quantum log-likelihood ratio of the model pe, it will be useful in a wide range of statistical estimation 
problems. 
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4.1 Achievability of the Holevo bound 

Theorem 12. Let {pe ; G 6 C M*^} he a quantum statistical model on a finite dimensional Hilbert 
space %, and fix a point 6q Cz 0. Suppose that pg ^ pg^ for all 6 £ 0, and that the quantum log- 
likelihood ratio Ch '■= £ (/3eo+/i|peo) is differ entiahle in h around h ~ and twice differ entiahle at 
h = 0. For any countable dense subset D o/R and any weight matrix G, there exist a sequence 
of estimators on the model ^ € that enjoys 

lim E^''^ [M(")] = h 

and 

lim TrG14^"^[M("'] - Cg, {pg,G) 

n— ^C30 

for every h Cz D. Here Cg„ {pg, G) is the Holevo bound at 9q. 

Proof. Let D :~ T^peg be the commutation operator with respect to the state pg^ (see Appendix]^, and 
let T be the minimal D invariant extension of the SLD tangent space spanjg {Li}"^^^ of the model {pg} 
■dt 6 ^ 9q, i.e., the smallest V invariant real linear subspace of Hermitian operators on Ji containing 
all the SLDs {iij^^i of pg at 9q. The minimality ensures that Tr pg^A = for all A e T because 
T' = {A e T ; Tr pggA ~ 0} is also 2? invariant. 

Let {Dj}^j^i be a basis of T, thus d < r. Let S be an r x r matrix whose (i, j)th entry is given 
by Yiij = Tr pggDjDi, and let r be an r x d matrix whose (i,j)th entry is given by = Tr pg^LjOi. 
According to Theorem 1191 in Appendix \X[ the Holevo bound for a weight G > can be expressed as 

Cgg{pg,G) = min{TrGZ + Tr VClmZVC ; Z ^ *FJ:F, 

F 

F is an r X d real matrix satisfying *FKe (t) = /}. (9) 

Letting 

1 " 

Xf") := ^ ^ /^^'^"'^ ® A ® (1 < i < r), 

k=l 

Corollary im asserts that ({pf "} , X^")) is jointly QLAN at 6*0, and that 

{x(-\pf;_^,^^)^N{{Rer)h,j:). (10) 

Let F be the matrix that attains the minimum in and let Z := *FS]F, V :~ KeZ, S := ImZ, 



V = \fG^ 



y/GlmZVG V(F^, and Z = V-y/^S. It then follows from CoroUarylMland Theorem 
]in Appendix [B] that 

Ggg ipg,G) = TTG(v + v) . 



Further, Lemma [T3l below assures that there exist a finite dimensional Hilbert space T-L and a state a 
and observables Bi {1 < i < d) on H such that Tr a Bi = and Ti a BjBi = Zij. Let 



where 1^") := F^X^'^ {l<i<d), 



4n) ._ 



In 



n 
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and / is the identity on %. A crucial observation is that (^x''"'\ pjf''^ , where p)^ 
converges to a classical Gaussian state: 



(") (gir. 



(x("\pi"))^7v(/.,v^+n 



(11) 



for all heW. In fact, 



lim Trpi"Mne^^'''- 



lim Tr p'"'^ 



= lim 



'Ph 



Tr/ 



Trcr 



v,t=l 



(12) 



\t=i 



\t=i 



where Xi := F^Xk (1 < « < d) are canonical observables with Xi, . . . , Xr being the basic canoni- 
cal observables of CCR(ImE) and {X,(f)h) ~ A^((Re r)/i, E), and Yi, . . . , are the basic canonical 
observables of CCR^^ImZ^ with {Y,i}j) ^ N{Q,Z). In the last line in we used as well as 

the quantum central limit theorem for F^"'. By using the explicit form ^ of the quasi-characteristic 
function for the quantum Gaussian state, (|12p is rewritten as 



\ r r \ r r 

t=l ^ ^ t=l s=t-|-l / \ " t=l t=l s=t+l 

= exp [j2 (v^Qh, - Uietiv + n^) -EE ^i^^iv + n 



t=l s=t+l 



This proves pi|) . 

Now according to Lemma [T4l below, there exist a quintuple sequence 



On 



of POVMs on [n®nj , taking values in a certain finite subset of M'', that enjoys the 

properties 

lim lim lim lim lim ;Ei"^[M("''"'^''?'P)] = /i, 



and 



p—^oo q—^oo foo m— ^C30 n— ^oo 



lim lim lim lim lim ' [M("^™'^'«'P)] F + 

p— fcxD g— >-oo ^— foo m— >-cx3 ri— ^oo 



r 



F(")r 



for all h e M'', where and denote the expectation and the covariance with respect to 

pI"^ . It then follows from Lemma [15] below that for any countable dense subset D of R'^ and any 
h Cz D, there exist a subsequence {{n,m{n),£{n)^q{n),p{n)}^^f^ such that 



lim = h, 



and 



This implies that the POVM Af ^"^ on "H®" that is uniquely defined by the requirement 



Trp^"^Af;^"^ = Tr ^p'"' (g) CT*^"^ ^^^n,m(n),£(n),g(n),p(n)) 
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for all density operator on -H®" and cj e f7(">™(")>'("),p(n),9(")) enjoys 

lim £;}"^[M(")] = h, 

n^oo 

lim K^"'[M(")] = y + f. 

n— f oo 

for all h E D. Recalling that Tr G{V + V^) = Cbq {pe, G), the proof is complete. □ 

Lemma 13. Given a d x d positive semidefinite Hermitian matrix J , there exist a finite dimensional 
Hilhert space % and a pure state p and ohservahles Ai (\ < i < d) on T-L such that Tt pAi = and 
Tr pAj Ai ~ Jij . 

Proof. Let U = C^+i, and let be a CONS of "H. We set |0) and Ylt^i [^j 1^) 

for i = 1, . . . , d. Then p := IV-) (V'l and A, (V-l + IV") {^ satisfy TrpA, = and pAjA, = 

Jij. □ 



Lemma 14. Given a sequence H*-"-* of finite dimensional Hilbert spaces, letX("^ = . . . , X^") 

6e a list of observables on "H*-"-*, and let family of density operators on "H^"-* parametrized 

byheWK If there is a real d x d positive definite matrix V such that 

^iV(/.,T/) (13) 

holds for all h £ W^, then there exist a quintuple sequence {M("^"'^'«'P);(n,m,^,5,p) G N^} of POVMs 
on that enjoy the properties 

lim lim lim lim lim [M("''"'^'«'P)] = /i, 



and 

p— >oo q~^oo >-oo ?n— >cxD ri— >oo 

Proof. Let 



lim lim lim lim lim F,^"^ [M("'™'^'«'P'] = K 

p— >cx3 q—^oo j-oo ?n— >cxD ri— >-oo 



?7i 2m J 

be a finite subset of R'', comprising {2m)'^ lattice points in the hypercube [—1,1]'^, and let f2('"'^'P) := 
p [-.p^p]d and ^2^™'^'^' := 17(™'«'P) u {0}. We introduce a Gaussian density function fL''\x) on 
M'' centered at a; = (wi, . . . , w^) S M'' by 

where x = (xi, . . . , Xd) G K'' and 

ai'Ht) (^) ' exp - , (s, t € M). 



By using this function, we define a POVM A/(",™^^9,p) = ^ g on "H*") that 

/(") 



takes values in the finite subset ft'^'^'^^ by 
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for w G and 



M, 



{n.m,£,q,p) 



(2to)' 



^(m,f,9)j-j^(ri)-j 



Here 



is the normalization with 



5W 



\/TTT 

Intuitively speaking, the dilFerence set fi'™^^' \ 17^'"'^'^^ works as a "buffer" zone that gives the default 
outcome ui = 0. This device is meaningful only when p < £. 
We shall prove that 



lim lim lim lim lim V P(w)Tr p|,"^Mi"'"''''«'P) = / P{uj)ph{uj)du 



(14) 



where P{lj) is an arbitrary polynomial of a; such that P(0) = and Ph{^) is a probability density 
function of the classical normal distribution N{h, V). Once has been proved, we can verify 

lim lim lim lim lim [M("''"'^'«'P)] = /i 



p— >oo g— >oo £— >oo m— >-oo ?i— >oo 



and 



lim lim lim lim lim Vr'[M^ 

p— ^oo q-^oo i—^oo m— ^C30 n— ^oo 



just by letting P{uj) = uji or P(a;) = UiUj (1 < i,j < d) in p^ . 
The first limit 7i oo in yields 

lim y P(w)Trpi"^Mi"^'"'^'«'P) 
= lim y P(cj)Trpi"^A4"^'^'^'«'P) 

n— foo — ^ 

^gO(™,f,p) 

= lim y P(tj)Trpi"^i?(™^^^9'(X(")) 

n— >cxD ^- — ^ 



/(") 



Ph{x)dx. 



(15) 



In the fourth line, we used the assumption and Corollary [51] in Appendix |Dl as well as the 
fact that functions gi'^\t) on K and g{t) on t > arc both bounded and continuous. Further, X = 
(Xi, Xd) is a classical random vector that follow the normal distribution N{h,V), and Eh[-] 
denotes the expectation with respect to N{h, V). As for the second limit m oo, due to 



P 



< max |P(w)| < oo. 
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the bounded convergence theorem yields 



Jim ([15]) = / hm T~(<i) 1 — Ph\x)ax 



where p^'^\uj, x) — (^) ^ cxp (^"f ~ "^i)^) : ^^id Darboux's theorem for the Riemann integral 

was used in the second line. Finally, the dominated convergence theorem and Fubini's theorem yield 

lim lim lim (fTB|) = lini lim / ^l. — j-j- pii{x)dx 



p->oo g-i-oo £-i.oo p-i-oo q^oo Jj^d Jj^a P^'^H^ , x) did 



p—¥00 q— >-00 



lim lim / ( / P{Lij)p''''\u!, x)duj \ ph{x)dx 



= lim lim / ( / p^*^ (w, x)p/i(x)dx ) P(aj)dw 

= lim lim / p^i^\LL>)P{uj)du! 
= lim / pii{uj)P{uj)duj 



p— f oo 



p,,(w)P(w)dcj, (17) 

where p^i^\uj) is the density function of N{h, V + ^I). This completes the proof. □ 

Lemma 15. For each i e N, let {a5j^„2---rvn' ("-I7 "2, . . . ,nr,n) e N(''+i)} be an {r + l)-tuple sequence 
on a normed space V . If, for each i € N, there exists an a* 6 V such that 

lim lim • • • lim lim a' = a\ 

ni— ^00 712—^00 n^.— ^00 n— >oo 

then there exist a subsequence {{ni{n),n2{n), . . . ,nr{n),n)}n^m that satisfies 



for all i G N. 



Proof. We first prove the case when r ^ 1. Let aj,^ := lim„_j.oo We construct a subsequence 

{(ni(fc), n(fc))}fceN in a recursive manner as follows. We set ni(l) = n(l) = 1. For fc > 2, it follows 
from lim„j^oo al^^ = q* that there exist an Ni{k) £ N such that ni > Ni{k) implies 

max a„, — a < — . 

l<i<k ' ^ ' k 

Thus the number ni(fc) max{A^i(fc), — 1) + 1} satisfies 



max 

l<i<k 



< 1^ (18) 



For this ni{k), it follows from lim„^oo a,\j(fc)„ = that there exist an N{k) £ N such that 

n > A''(fc) implies 

1 



max 

i<i<fc 



< -. (19) 
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Thus we set n{k) := max {7V(fc), ri(fc - 1) + 1}. 

Now let k(n) :~ max{fc; n{k) < n}, which is non-decreasing in n and hm,i_5.oo k{n) = oo. We show 
that the subsequence {ni{k{n)),n); n e N} enjoys the required property: for all i e N, 



lira a] 

n— f oo 



ni{k{n))n — ^ 



Given i e N and e > arbitrarily, there exist an iV e N such that n> N implies k{n) > max { j, [|] }. 
Then for all 7i> N, we have 



< 



^ni(fc(n))n '^ni(fc(n)) 



"ni(fc(n)) ~ 



3 J 
ni{k{n))n ni{k{n)) 



+ max 
i<j<fc(«) 



ni(fc(n)) 



< max 

l<j<k{n) 

In the third inequality, we used (|T8|) and (fT9|) . as well as its premise n > n{k{n)) > N{k{n)). 
The proof for a generic r is similar. 



□ 



4.2 Application to qubit state estimation 

In order to demonstrate the power of our method, we explore qubit state estimation problems. 
Example 16 (3-dimensional faithful state model). 

The first example is an ordinary one, comprising the totality of faithful quit states: 



S{C- 



e e 



where cr^ (i ~ 1,2,3) arc the standard Pauli matrices and 8 is the open unit ball in M.^. Due to the 
rotational symmetry, we take the reference point to be 9q = (0,0, r), with < ^ < 1- By a direct 
calculation, we see that the SLDs of the model pe &i = Oq are (Li, L2, -^3) = (ci, (J2, {rl + 0-3)"^); 
and the SLD Fisher information matrix J^"^' at 6q is given by the real part of the matrix 




Given a 3 x 3 real positive definite matrix G, the minimal value of the weighted covariances at 
9 ^ 9o is given by 

minTrGl/e„[M] = {pe,G), 

M 

where the minimum is taken over all estimators M that are locally unbiased at 6q, and 

2 

-(1) 



is the Hayashi-Gill-Massar bound [II [3] (see also [21]). On the other hand, the SLD tangent space is 
obviously 2? invariant, and the Holevo bound is given by 

Cg„ {pe, G) := Tr GJ^"^'' + Tr ^/G Im ,/(-^)"' ^/G 
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where 



(Rc J)-V(Rc J)-^ 




is the nwerse RLD Fisher information matrix (See Coroharv [20l in Appendix 

It can be shown that the Hayashi-Gill-Massar bound is greater than the Holevo bound: 

Let us check this fact for the special case when G = j'"^'. A direct computation shows that 



and 



= 3 + 2r. 



The left panel of Figure [T] shows the behavior of Ceg [pe, J'"^^) (sohd) and Cg^J [pe, J^"^') (dashed) as 
functions of r. We see that the Holevo bound Ce^ {pe, J^^^) is much smaller than Cg^J (pe, J^"^)). 

Does this fact imply that the Holevo bound is of no use? The answer is contrary, as Theorem [121 
asserts. We will demonstrate the asymptotic achievability of the Holevo bound. Let 



1 " 



k=l 



and let X,^"'' :— A^-"'' for i = 1, 2, 3. It follows from the quantum central limit theorem that 



^ N \0, 



J J' 
J J, 



Since 



m ■.= Cipe\peo)^2\og 



is obviously of class C°° in 0, Corollary [TT] shows that ({pf "} , X'")) is jointly QLAN at 6*0, and that 

for all /i G M'^. This implies that a sequence of models ^pf"^^/^ i ^ <= converges to a quantum 
Gaussian shift model {iV((Re J)/i, J) ; /i e M'^}. Note that the imaginary part 



S 



( ° 
V 





oy 



of the matrix J determines the CCR(S'), as well as the corresponding basic canonical obscrvablcs 
X ~ {X^,X'^,X^). When r 7^ 0, the above S has the following physical interpretation: and 
X^ form a canonical pair of quantum Gaussian observables, while X^ is a classical Gaussian random 
variable. In this way, the matrix J automatically tells us the structure of the limiting quantum 
Gaussian shift model. 
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Figure 1: The left panel displays the Holevo bound C(o,o,r) (pe,"/^'^') (solid) and the Hayashi-Gill- 
e, J^^^) (dashed) for the 3-D model pe ^ ^ {l + O^ui + 6*20-2 + O^a^) as func- 



Massar bound C, 
tions of r = 



(1) 
(0,0,1 



bound C, 



(1) 

(0,r) 



The right panel displays the Holevo bound C(o.r) (pe, J^'^') (solid) and the Nagaoka 
(pe, J^^^) (dashed) for the 2-D model pg = 



Now, the best strategy for estimating the shift parameter h of the quantum Gaussian shift model 
{A^((Re J)h, J);h&W^]\s the one that achieves the Holevo bound Ch (iV((Re J) h, J),G), (of., The- 
orem[2n]in Appendix IbI). Moreover, this Holevo bound Ch (iV((Re J) ft., J), G) is identical to the Holevo 
bound C^o (pe, G) for the model pg at 6*0, (cf.. Corollary [531 Recall that the matrix J is evaluated 
at of the model pg). Theorem [T2l combines these facts, and concludes that there exist a sequence 
M*^"^ of estimators on the model h € M'^j that is asymptotically unbiased and achieves 

the common values of the Holevo bound: 

lim TrGV^^^"^[Af = Ch {N{{Re J)h, J),G) = Cg„ {pg, G) 

for all h that belong to a countable dense subset of R'^. 

It should be emphasized that the matrix J becomes the identity at the origin 6*0 = (0,0,0). This 
means that the limiting Gaussian shift model {N{h, J) ; /i e M^} is "classical." Since such a degenerate 
case cannot be treated in [51[TTJ|S], our method has a clear advantage in applications. 

Example 17 (Pure state model). 

The second example is to demonstrate that our formulation allows us to treat pure state models. 
Let us consider the model S = {|V'(6'))(^(6')| ; = (6l*)i<,<2 £ &} defined by 

1 i(eVi+9^(T2) 




Y^cosh \\9" 

where Q is an open subset of containing the origin, and j| • j| denotes the Euclid norm. By a direct 
computation, the SLDs at 9q = (0, 0) are (ii, L2) = {ci, (T2), and the SLD Fisher information matrix 
j'"^) is the real part of the matrix 



J = [Tr 



that is, J^"^) — I. Since the SLD tangent space is V invariant [T], the Holevo bound for a weight G > 
is represented as 

Cg„{pg,G) :=TrGj(^)"' + Tr VClm J^^'>'' VG 
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where 

( 1 -V^^ 



-1 1 



is the inverse RLD Fisher information matrix (see Corollary in Appendix IX)) . 
Let us demonstrate that our QLAN is applicable also to pure state models. Let 



1 " 

^ k=l 



and let := A^-"'' for i = 1,2. It follows from the quantum central limit theorem that 



iV 0, 



J J , 



Since 

£(6) -.^ Cipg\pg,) = e^ai + e^a2 - \ogcosh\\e\\ 

is of class C°° with respect to 0, it follows from Corollary [IT] that {{pf"} , is jointly QLAN at 

00, and that 

(^^"\C+VV^) ^ Ar((Re J)/., J) = J(«)-') 
for all /i e M^. Theorem [T2l further asserts that there exist a sequence M*^") of estimators on the model 
\^pf\h/^ ; /i G that is asymptotically unbiased and achieves the Holevo bound: 

Jim^TrGF^ = C,, (7V(/i, g) C(o,o) (pe, G) 

for all h that belong to a dense subset of M.^. In fact, the sequence can be taken to be a separable 
one, making no use of quantum correlations [15]. (See also Appendix [C] for a simple proof.) Note that 
the matrix j'^^'^ ^ is degenerate, and the derived quantum Gaussian shift model ^N{h, J^^'f ^)| is a 
canonical coherent model [T]. 

Example 18 (2-dimensional faithful state model). 

The third example treats the case when the SLD tangent space is not 2? invariant. Let us consider 
the model 

5 = |p0 = i (/ + eVi + + zoVl- II^IP^s) ; 9 = {0')i<:<2 e e| , 

where < zo < 1, and is the open unit disk. Due to the rotational symmetry around z-axis, we 
take the reference point to be 6*0 = (0,r-), with < < 1- By a direct calculation, wc see that the 
SLDs at 0Q are (ii, L2) = ^ci, ^^{(72 — rl)j . It is important to notice that the SLD tangent space 
span {Li},^^j^ is not T) invariant unless r = 0. In fact 

Vai = z{r)a2 - ra^, T)a2 ^ ~z{r)ai. 



where z{r) := E[a3] = zq\J\ — . The minimal T) invariant extension T of the SLD tangent space has 
a basis D2, D^) :~ (Li, L2, (T3 — z{r)I). The matrices E, J. and r appeared in Definition [5] and 
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Corollary [TT] are calculated as 



1 



z{r) 
T- z(r) 



z(r)2 



z{r) 



z{r) 



1 



I 



1 



z{r) 



V z{r) 
1 



0(r)2 / 



[Trpe„LjCr,].^. 



z{r) 



-rv/— 1 — z(r 



z{rY 



r 
z{r) 



1 zt 



Given a 2 x 2 real positive definite matrix G, the minimal value of the weighted covariances at 
6 ^ 6q is given by 

minTrCyeolM] -d'^ {pe,G), 

M 

where the minimum is taken over all estimators M that are locally unbiased at 9q, and 

2 

.(1), 



is the Nagaoka bound [T6| (see also j2T]). 

It can be shown that the Nagaoka bound is greater than the Holevo bound: 

d,'J ipg,G)>CeAPo,G). 

Let us check this fact for the special case when G = j'"^'. A direct computation shows that 



C, 



and 



G, 



2(l + zo)-r2(l~z2), iiO<r< 

V i — 



if 



1 _ -.2 

1 Zq 



< r. 



The right panel of Figure [T] shows the behavior of Ceg (pe, ^Z'-'^-') (solid) and Gg^^ (pe, ^Z'-'^-') with 
zq ~ J (dashed) as functions of r. We see that Holevo bound Geo {potJ^^^) is much smaller than 

As in Example [T51 we demonstrate that the Holevo bound is asymptotically achievable. Let 

a|") := ^ y I"""-' ®L,® {i = 1, 2), 
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and let 

xf) :=^y/«'=-l®z?,®/«"-^ (j = 1,2,3). 

It then follows from the quantum central limit theorem that 

(C:;:.Vr),«Hf. ;)). 

Therefore, Corohary [lU shows that ({p^"} , is jointly QLAN at 6*0, and that 

for all h(^M?. 

It should be noted that the off-diagonal block t of the "quantum covariance" matrix is not a square 
matrix. This means that the derived quantum Gaussian shift model {Af((Re t)/i, S) ; h e M^} forms 
a submanifold of the total quantum Gaussian shift model derived in Example 1161 corresponding to 
a 2-dimensional linear subspace in the shift parameter space. Nevertheless, Theorem [T^] asserts that 
there exist a sequence M("' of estimators on the model {p^^," ^/y^^; /i G K'"^! that is asymptotically 
unbiased and achieves the Holevo bound: 

lim TrGl/,["'[M(")] = Ch (iV((Re r)/i, S), G) = Ce„ (pe, G) 

n— f oo 

for all h that belong to a dense subset of R'^. 

4.3 Translating estimation of h to estimation of 6 

As we have seen in the previous subsections, our theory enables us to construct asymptotically optimal 
estimators of h in the local models indexed by the parameter 9q + h/^/n. In practice of course, Oq is 
unknown and hence estimation of /i, with 9q known, is irrelevant. The actual sequence of measurements 
which we have constructed depends in all interesting cases on Oq. 

However, the results immediately inspire two-step (or adaptive) procedures, in which we first mea- 
sure a small proportion of the quantum systems, in number ni say, using some standard measure- 
ment scheme, for instance separate particle quantum tomography. From these measurement outcomes 
we construct an initial estimate of 6*, let us call it 0. We can now use our theory to compute the 
asymptotically optimal measurement scheme which corresponds to the situation 6q — 9. We proceed 
to implement this measurement on the remaining quantum systems collectively, estimating h in the 
model 9 = 9 + h/y/n2 where n2 is the number of systems still available for the second stage. 

What can we say about such a procedure? Ifni/n— >-a>Oasn— >-oo then we can expect that the 
initial estimate 9 is root n consistent. In smooth models, one would expect that in this case the final 
estimate 9 = 9 + hj ^fn^ would be asymptotically optimal wp to a factor 1 — a: its limiting variance 
will be a factor (1 — a)^^ too large. 

If however ni — ^ cxd but ni/71 — > a = then one would expect this procedure to break down, 
unless the rate of growth of ni is very carefully chosen (and fast enough). On the other hand, instead 
of a direct two-step procedure, with the final estimate computed aa 9 + hj ^pwl. one could be more 
careful in how the data obtained from the second stage measurement is used. Given the second step 
measurement, which results in an observed value /i, one could write down the likelihood for h based on 
the given measurement and the initially specified model, and compute instead of the just mentioned 
one-step iterate, the actual maximum likelihood estimator of 9 based on the second stage data. Such 
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procedures have earlier been studied by Gill and Massar [5] and others, and shown in special cases to 
perform very well. 

However, in general, the computational problem of even calculating the likelihood given data, 
measurement, and model, is challenging, due to the huge size of the Hilbert space of n copies of a finite 
dimensional quantum system. 

5 Concluding remarks 

We have developed a new theory of local asymptotic normality in a quantum regime based on a 
quantum extension of the log-likelihood ratio. This formulation is applicable to any model satisfying 
a mild smoothness condition, and is free from artificial setups such as the use of a special coordinate 
system and/or non-degeneracy of eigenvalues of the reference state. We also have proved asymptotic 
achievability of the Holcvo bound for the local shift parameter on a dense subset of the parameter 
space. 

There are of course many open questions left. Among others, it is not clear whether every sequence 
of statistics on a QLAN model can be realized on the limiting quantum Gaussian shift model. In 
classical statistics, such a problem has been solved affirmatively as the representation theorem, which 
asserts that, given a weakly convergent sequence T'"^ of statistics on ; h G there exist 

a limiting statistics T on {N{h, J^^) ; h G R"*} such that T^"^ ^ T. Representation theorem is useful 
in proving, for example, the non-existence of an asymptotically superefficient estimator (the converse 
part, as stated in Introduction). Moreover, the so-called convolution theorem and local asymptotic 
minimax theorem, which are the standard tools in discussing asymptotic lower bounds for estimation 
in LAN models, immediately follows |20| . Extending the representation theorem, convolution theorem, 
and local asymptotic minimax theorem to a quantum regime is an intriguing open problem. However 
it surely is possible to make some progress in this direction, as for instance the results of Gill and Gu^a 
g]. In that paper, the van Trees inequality was used to derive some results in a "poor man's" version 
of QLAN theory; see also [2]. 

It also remains to be seen whether our asymptotically optimal statistical procedures for the local 
model with local parameter h can be translated into useful statistical procedures for the real world 
case in which 9q is unknown. 

Appendices 

A Commutation operator and the Holevo bound 

In the study of quantum statistics, Holevo [9] introduced useful mathematical tools called the square 
summable operators and the commutation operators associated with quantum states. Let "H be a 
separable Hilbert space and let p be a density operator. We define a real Hilbert space C\{p) associated 
with p by the completion of the set BhiT-L) of bounded Hermitian operators with respect to the pre- inner 
product (AT, Y) ^ ReTrpA'y. Letting p ~ Sj\il>j){il}j \ be the spectral representation, an element 
X e C\{p) can be regarded as an equivalence class of those Hermitian operators, called the square 
summable operators, which satisfy '^ill-'^V'jP < oo (so that -^j S Dom(X) if Sj ^ 0) under the 
identification Xi ~ X2 if Xiijjj = X2ipj for Sj ^ 0. The space C\{p) thus provides a convenient tool to 
cope with unbounded observables. Note that when % is finite dimensional, the setup is considerably 
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simplified to be Cf^{p) = BhiV.) / ker p. 

Let C^{p) be the complexification of C\{p), which is also regarded as the completion of B{'H) with 
respect to the pre-inner product 

\tvp{YX* +X*Y). 

Thus C^{p) is a complex Hilbert space with this inner product. Let us further introduce two sesquilinear 
forms on BCH) by 

{X,Y)p -.^Ti-pYX*, [X,Y]p ^^Ty p{YX* - X*Y). 

and extend them to C^{p) by continuity. 

The commutation operator Dp : C^{p) — > €?{p) with respect to p is defined by 

[X,Y]p^{X,VpY)^, 

which is formally represented by the operator equation 

Vp{X)p + pVpiX) = ^(Xp- pX). 

(To be precise, Holevo's original definition is different from the above one by a factor of 2.) The operator 
"Dp is a C-linear bounded skew-adjoint operator. Moreover, since the forms [ • , • ]p and ( • , • )p are real on 
the real subspace C\{p), this subspace is invariant under the operation of Dp. Thus Dp can be regarded 
as an M- linear bounded skew-adjoint operator when restricted to C\{p) as Vp : C\{p) — > C\{p). When 
no confusion is likely to arise, we drop the subscript p of Vp and simply denote it as P. 

Let S = {pe; 9 G Q G K''} be a quantum statistical model satisfying the conditions: 1) the 
parametrization 9 i-^ pg is smooth and nondcgcncratc so that the derivatives {dpg/d9'}i<ci<d exist in 
trace class and form a linearly independent set at each point 9 G O, and 2) there exists a constant c 
such that 



d 



<c{X,X) 



Pe 



for all X e B{7i) and i. The second condition assures that the linear functionals X i-> {d/d9'^)Tr pgX 
can be extended to continuous linear functionals on C^(pg). Given a quantum statistical model sat- 
isfying the above conditions, the symmetric logarithmic derivative (SLD) Lg^i in the ith direction is 
defined as the operator in {pp^^ ) satisfying 



-^TxpgX = {Lg^,,X)p^. 



It is easily verified that Lg^i e L\{pg)] so the definition is formally written as 

1 



^{Lg^iPg + pgLg^i). (20) 

When no confusion occurs, we simply denote Lg,?: as Li. Since Li is a faithful operator representation 
of the tangent vector 5/96'*, we shall call the M-linear space span^{Li}f^j^ the SLD tangent space of 
the model pg at 9. Incidentally the d x d real symmetric matrix Jg := [KeTr pgLiLj]-^^^- is called 
the SLD Fisher information matrix of the model S at 9. 

An estimator M for the parameter 9 of the model S is called unbiased if 

Eg[M]=9 (21) 
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for all 6* G 8, where Eg[-] denotes the expectation with respect to pg. An estimator M is called 
locally unbiased at 9q Cz O ii the condition (j2ip is satisfied around 9 = Oq uj> to the first order of the 
Taylor expansion. It is well known that an estimator M that is locally unbiased at 6*0 satisfies the 
quantum (SLD) Cramer- Rao inequality, Vg^ [M] > Jg^ , where Vg^ [ ■ ] denotes the covariance matrix 
with respect to pe^ . The lower bound Jg^ cannot be attained in general due to the non-commutativity 
of the SLDs. Because of this fact, we often switch the problem to minimizing the weighted sum of 
covariances, TrGV6/o[M], given a d x d real positive definite matrix G. It is known that this quantity 
also has a variety of Cramer- Rao type lower bounds [S] : 

TTGVg„[M]>C9„ {pg,G). 

Among others, we concentrate our attention to the Holevo bound [S]: 

Cg„ {pg, G) min{Tr GV ] V is a. real matrix such that V > Z{B), Z,j (B) ^ Tr pg^BjE^, 

Bi, . . . , Bd are Hermitian operators on H such that RcTt pggLiBj ~ Sij}. (22) 

The minimization problem over V is explicitly solved, to obtain 

Cg,{pg,G) = mm{Ti-GZ{B)+Ti- VGlmZ{B)VG ; Z^j{B)^Ti-pe„B,B„ 

B 

Bi, . . . , Bd are Hermitian operators on % such that ReTv pg^LiBj ~ 5ij}. 

Our aim here is to derive a further concise expression for it in terms of the minimal 2? invariant 
extension of the SLD tangent space. 

Theorem 19. Given a quantum statistical model {pe ; G O C M'*} on H, let T be the minimal T) 
invariant extension of the SLD tangent space spanjj{Li},f^j^ of the model at 6 ^ Oq, and let {Dj}^^i 
be a basis ofT- The Holevo bound defined by iSS^) is rewritten as 



Gg„ ipg,G) = min{ TtGZ + Tr 



GIuiZVG 

F 

t 



FEF, 



F is an r X d real matrix satisfying FUe (r) = /}, (23) 

where E and r are rxr and rxd complex matrices whose {i,j)th entries are given by = Tr pg^DjDi 
and Tij = Tr pgf^LjDi. 
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Proof. Let be the orthogonal complement of 7~ in £^{pgg) with respect to the inner product 
and let V : C\{pg„) — > T and V-^ : C\{pgg) — > be the projections associated with the decomposition 
Cl{pg„) = r ® T-L. Note that if X e and Y eT, then 

We show that the operators {Bj}'^^^ that achieve the minimum in can be taken from T. Let 
{Bj}'^^^ C £^(peo) satisfies the local unbiasedness condition KcTt pg^LiBj — Sij, which is rewritten 
as 

Then {V{Bj)}j^^ also satisfies the local unbiasedness 

iL.,V{B,))^{L,,B,) =S.,. 
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Further, 



- {V{B.,IV{B,))^^^ + {V^{B,),V^{B,))^^^^ ^ Z.,,{V{B)) + Z.,,{V^{B)). 

Since Z( • ) is a Gram matrix and is positive semidefinite, this decomposition imphes that Z{B) > 
Z(V{B)). Thus the obscrvables B that minimize can be taken from T. 

Let Bj € T be expanded as Bj ~ FjDk, where F is an r x d real matrix. Then the local 
unbiasedness condition is rewritten as 

or in a matrix form, 

*F(ReT) =/. 
Further, the Gram matrix Z{B) is rewritten as 

Z,,iB) = {B,,B,)^^^ = Fj'F^HD,,D,)^^^ , 

or, 

Z{B) = *FSF. 

This proves the claim. □ 

When the SLD tangent space itself is D invariant, the Holevo bound can be represented in terms 
of the RLD Fisher information matrix as follows. 

Corollary 20. Let {pe ; 6* e G C M^} be a quantum statistical model, and let Li (1 < i < d) be the 
SLDs at Oq. If the SLD tangent space spang {Li}'^^^ at Oq is D invariant, then 

Ceo (P8,G) =TrG(j(^')-i +Tr VGIm{j'^'^^)-^VG 
where := (Re J)"^ J (Re J)~^ with = Ti pe„LjL,. 

Proof. Let us set Di := Li for 1 < i < d in Theorem [T9l Then S = t, and the local unbiasedness 
condition *F (Rer) = / has a unique solution F = (Re E)"', whereby Z = (Re J) ,J (Re J) . □ 

Note that RLDs may not exist if the model is degenerate (i.e., non- faithful) . This means that 
may not be well-defined for such a model. Nevertheless we use the notation (j'^^) ^ even for a 
degenerate model, and call it the inverse of the RLD Fisher information matrix, as long as the SLD 
tangent space is T) invariant. For an idea behind this nomenclature, consult [T]. 

Finally, we show that the Holevo bound for the nth i.i.d. extension model is precisely i times that 
for the base model. 

Corollary 21. Given a quantum statistical model S ^ {pe ; 6* e 9 C W^], letS^'''^ = {pf" ; 6* e 9 C M'*} 
be the nth i.i.d. extension model. Then 

Ceo {pT\G) = -Ceo ipe,G) . 
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Proof. Let us distinguish quantities that belong to models of different extension by specifying the 
degree n of extension in the superscript. Letting {Li}f^^ and {Dj}^^^ be SLDs and a basis of T in 
Theorem [191 the corresponding quantities for 5^"^ are given by 



and 



Thus 



so that 



and 



fc=i 



doe to Theorem [TH □ 

B Estimation of quantum Gaussian shift model 

In this section, we briefly overview the estimation theory for a quantum Gaussian shift model. For a 
mathematically rigorous treatment, consult |^. 

Lemma 22. Let (X, (ph) ^ N{h, J), where J is a d x d positive semidefinite complex matrix. Then 

MX^)^h, (24) 

and 

4>h{{Xj-h,){X,~hi))^J,, (25) 

hold. 

Proof. Letting U{£,) := e^«"^', 

Mum = i + v^Mcx^)-lM{ex,f) + o{e) 
= 1 + v^Mx^)e - lMx^x,)ce + o{e) 
= 1 + v^Mx^)c - \Mx^ o x/jce + o{e), 

where oXj^\ {XtXj + XjXi). Further, letting F = Rc J and S* = Im J, 

1 ... ... 



= 1 + v-if/i. - (v^y- + h.^i)i'i' + o(r ). 



2 
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A comparison immediately leads to and the identity (j>hiXi o Xj) ~ Vij + hihj. Thus 

MiXj ' h,)ix, - h,)) = 4)H{x,x,-hjX,-h,Xj + h,h,) 

= (t>h{XjXi) ~ h^hj 



~H X, o Xj - -[X,, X,] - h,h, = X, 



□ 



In what follows, we treat the quantum Gaussian shift model {iV(r/i, S) ; /t. G M''} on CCR (ImE), 
where E is an r x r complex matrix such that S > and Re E > 0, and r is an r x c? real matrix with 
d<r such that rankr = d. Let X = {Xi, . . . , Xr) be the basic canonical observables of CCR(ImE), 
and {X,(l)h) ^ N{Th,T,). 

Lemma 23. Let U{^) c^^'^- . The SLD L, (1 < i < d) at h defined by 



is given by 



r 

Lfe = V[(ReE)-Vl {X,~{Th),I). 



(26) 



(27) 



Proof. In this proof we lift Einstein's summation convention. Let V = ReE and S = ImE, and fix 
fc e {1, . . . , d} arbitrarily. Due to the Baker- Hausdorff formula, 



U{0 = e^SZLiTx, = exp (^"^E ^fe^^'^'j ^xp cxp 
By differentiating in .f'"', we have 

^^u{0 = -^i(jZsuC-x^ uio. 



Thus 

MiXk - {Th)ki)u{0) 



hJlj2s,.e-v^^^irh),i\uio 



Vi=i 

^ r 



-1 



-1- 



d 

A- 



^ {Th)k C^ 



Similarly, we obtain 



MuiOiXk - {Th)ki)) = V^iJOkMum- 



(28) 
(29) 
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By combining ^ and (HH), 



( [Xk - iTh)ki)uiO + uiOiXk - {Th)ki) ) = 2V^{vOk Mum- (30) 

On the other hand, by a direct calculation 



Q^-MUiO) = g^c^-^ ^^"-^ = ^('eT)fc0.(C/(O). (31) 



A comparison between (|30|) and ([31]) yields 

r 

^fc = E[^"M^fc (X, ~ (r/i),/). 



f=i 



□ 



Let Lk := Xfe - {Th)kl. It follows from ^ and (EH) that V^^{L,) = EI=i(^"^'S')fe»-£fe, where 
is the commutation operator with respect to 4>h defined by 

{u{Ov^, {X) + v^, {x)u{0) = V^4>h {u{Ox - xu{0) ■ 

This means T = span ■{ ife f is I?^,^ invariant. Further, we can check from (j27p that span {iij^^j^ C 
T and 

4>h{LjLi) = (32) 

and 

Re0/i(ijXi) = n-j. (33) 

These relations play a fundamental role in connecting a general quantum statistical model S = 
{jOe ; S C M'*} on % with a quantum Gaussian shift model Q = {N{Th, Y,) ; h £ R''} as follows. 
Let {£f be the SLDs of the model 5 at 6* 6*0, and let the minimal invariant extension 
of the SLD tangent space span{Lf}f^^. Further let {D'j}^^^ be a basis of and let S and r are 
r X r and rxd matrices whose (i, j)th entries are given by = Tt pg^DjOi and = KcTi pg^LjOi. 
Based on those information, we introduce a quantum Gaussian shift model Q — {iV(r/i, E) ; h dW^^ 
on CCR(ImE), which exhibits relations ([5^ and ([55)1 . Recall that the Holevo bound of a quantum 
statistical model is completely determined by the information E and r (Theorem [T9| . We thus obtain 
the following important consequence. 

Corollary 24. The Holevo bound Cq^ {pg, G) for the model S at 9 — 9q is identical to the Holevo 
hound Ch {N{Th, E), G) for the Gaussian shift model Q . 

As to the achievability of the Holevo bound Ch {N{Th, E), G) for the Gaussian shift model G, we 
have the following. 

Theorem 25. Given a weight G > 0, there exist an unbiased estimator M that achieves the Holevo 
bound for the model {N{Th, E) ; /i G R"*}, i.e., 

TTGVh[M]^Ch{N{Th,J:),G). 

Proof. Let F be the matrix that achieve the minimum of for the model {A^(r/i, E)}/j, and let 
Z = ^FT.F. Further, let V = ReZ, S = lmZ. V = VcF^ VGIuiZVG VcF^ , and Z = V - ^f^S. 
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Wc introduce an ancillary quantum Gaussian state {Y, ip) ~ (0, Z) on another CCR , and a set 

of canonical observables 

'Xi:=Xi®I + I'^Yi {l<i<d), 

on CCR (g) CCR|^— S*^, where Xi ~ F^X^. It is important to notice that the CCR subalgebra 
A[X] generated by {Xi}i<i<d is a commutative one because 

[Xi,Xj\ = - Si J = 



for 1 < i, j < d. Moreover 



This means that the observables Xi {\ < i < d) follow the classical Gaussian distribution N{h, 1/ + F). 
In particular, 

= h 



for all /i e M'', and 



Tr GVh [X]=Tt:G{V + V) =Ch{N{Th,Y.),G). 



The claim was verified. 



□ 



C Estimation theory for pure state models 

Lemma 26. Let p he a pure state and Ai, . . . , Ad observables on a finite dimensional Hilbert space 
%. If Jij := TrpAjAi are all real for 1 < i, j < d, there exist observables Ki, ... , Kd such that 

[A,+ K,,Aj + Kj] = 0, 

for 1 < i, j < d and 
for 1 < i < d. 



Proof. Let p and let \li) := Ai for 1 < i < d. Because ("01/^) and {li\lj) (= J,-,;) are all 

real, there exi 
for all k. Let 



real, there exist a CONS {\ek)}'l™\^ of H such that (ehlip) and {e^lli) are all real, and that (efc|?/') 7^ 



dim y. / I , \ 

^' •= g 

and K, := Ai - Ai. Obviously [A^ + Ki, Aj + Kj] = [Aj, A^] = 0, and 

10) = (i, - A,) 10) = |/,) - |/,) = 0. 
This means Kip = 0. □ 

Theorem 27. Let {pe ', G C M'^} &e a quantum statistical model comprising pure states on a finite 
dimensional Hilbert space H, and let Cbq [pe, G) be the Holevo bound at9o G Q for a given weight G > 0. 
There exist a locally unbiased estimator AI at 9q G Q such that TrGy[A/] = C^q ipe,G). 
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Proof. Let T be the mmimal D invariant extension of the SLD tangent space span{Li}^^-^ of the 
model {pe} at 9 = 6q, i.e., containing all the SLDs {Li}'^^^ of {pe} at 6*0, let {DjYj^i be a basis of T. 
Let E, r be r X r, r X d complex matrices defined by Y^j = Tt pg^DjDi, Tij = Tv pg^^LjDi. According 
to Theorem [HI the Holevo bound for a weight G > can be expressed 



Coa [pe, G) 



minjTr GZ + Tr 

F 



GlmZVG 



Let F be the matrix that attains the minimum in (|34p . and let Z 



V 



GlniZVG 



V(f^, and Z = V 



; z = *fj:f, 

F is an r X d real matrix satisfying *FRe (r) = /}. (34) 

*FJ:F, V ReZ, S := ImZ, 
IS. Lemma [13] assures that there exist a Hilbert 



space H and a pure state cr and observables Bi (l < i < d) on H such that Tr aBi = and Tr aBjBi 



Zij. Further, let 



X, ■■= X,(g) i + Kg) B^ {l<i<d), 



where Xi := F^D^ {I < i < d), and / is the identity matrix on "H. It then follows that 



Tr ipg„®c7)X,X, = {V + V 



(35) 



According to Lemma [26l there exist observables Ki, . . . , Kd on T-L®Ti such that [Xi + Ki, Xj+Kj] = 



and Ki {pe^ ® a) —0. Let Ti 



\I®I+{X, + Ki). ThenTi, 



Td are simultaneously measurable. 



and satisfy the local unbiasedncss condition; 

Tr [pdo 

and 

Tr {d,pe, ® o) f, 



Further 

Ve„[f],,j = Tr [pg^ 
This completes the proof. 



= Tid^peoXj 

= F^^TTd,p9„Dk 

= Fj^ReTr pe,MDk 

= {F(Rer)},, = 5,, 



a) {X^ + K,) {X, + lu) = (V + V 



□ 



D Quantum central limit theorem 

Jaksic, Pautrat, and Fillet |10| proved the following strong version of a quantum central limit theorem. 

Proposition 28. Given a sequence -H^") of Hilbert space, let p^'^^ and A^") = (^i"\ • ■ • , ^i"^) « 
state and a list of observables on "H^"' that enjoy the quantum central limit theorem in the sense of 
convergence of the quasi- characteristic function: 

■^N{h,J)^{X,cj,), 

where J is a dx d positive semidefinite matrix. Then for any bounded continuous functions /i, . . . , fm 
and a noncommutative polynomial P, it follows that 

Jim Trp(")p (7^) =^{P (JV^)) . 
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where J(b] := ... , fi{Bd), . . . , /,„(Si), . . . , f^{Ba)) for a given list B = (Bi, . . . , Bd) of 

observables, and P (JiB)) := . . . , fi{Bd), . . . , fM)). 

Proposition[28]is strong enough to prove the foUowing, which is essential in constructing a sequence 
of POVMs that asymptotically achieves the Holevo bound (Theorem [T2|) . 

Corollary 29. Under the same assumption as in Provosition \ 28\. for any bounded continuous func- 
tions g,fi, . . . , fmj o,iT-d- noncommutative polynomials P,Q, with P being Hermitian operator-valued, 
it follows that 

Jim Trp("),g (^P (^./(A^) ) Q g {p 



g[P[^f{X)))Q[^f{X)) g[P[^f{X 

Proof. Let I := maxi<i<m sup^ There exist Ip > Q and /q > such that Ip > P(~^) 

and Iq > QC^) for any list ~^ ~ Bdm) of observables such that < /. Let Ig := 

sup{|g(a;)| ; x £ [-~lp,lp]}. There exist a sequence R^'^\x) of polynomials that uniformly converges 
to g{x) on [-lp,lp]. 
Let 

:= Trp(")i?W (^P (J(A^)^^ Q (j(A^)^ R^'''> (p (j{A^)^^ , 



a-kn 



and let 



a„:=TrpMg(^P i^fiA^"^) 
We show that akn uniformly converges to a„ as /e — )- oo. In fact, letting In :— sup {R'^'^^x) ; k £ N,x £ [—If 



sup |a„ - ttkn 
nGN 



sup 

raeN 



Trp(")g (^P (JiA^)^'^ Q (JiA^)^ g (^P (7{A^) 
Trp(")i?W (^P (Ji^)^^ Q (^f{A^)^ i?"^-) (^P (^f{A^)^ 



< sup 

rteN 



+ sup 

nGN 



< Iglq sup 



Trp("\9(p(/(A("')) ) Q(/(A("') 
Trp(") 



nGN 

-\-IqIp, sup 

nGN 



g (^P (J{A^)^ ^ - -R^^' (^P (J(A^)^ ^ Q {J{A^)^ R^^^ {p (^J{A^) 
g (^P - i?W (^P 
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Pif{A^"^)] ] -R'-''^ (p( /(A^")) 



<Iq{Is + Ir) sup g{x) ~ R'^''\x) 
xe[-ip,ip] 

which converges to zero as k oo. 

The uniform convergence akn =^ in as well as the existence of limfc_yoo hm„_j.oo flfcn, which follows 



33 



from Proposition [211 ensure that 

Jim Trp(").g (^P (^/(I^) ) Q (^J{A^)^ g (^P (^/(A^"'))) 

= lim Km Trp^'^^i^^'^) (^P (j(A^)^^ Q {j{A^)^ R^''^ (^P (7{A^)^^ 
= lim lim akn 

n— >oo fc— >oo 

= lim lim a^n 

k—^oo n— ^oo 



= (5 j Q [fix)) 9 [P [fix) 

This proves the claim. □ 
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